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Abstract 

In  this  work,  a  one-dimensional,  isothermal  two-phase  mass  transport  model  is  developed  to  investigate  the  water  transport  through  the  membrane 
electrode  assembly  (ME A)  for  liquid-feed  direct  methanol  fuel  cells  (DMFCs).  The  liquid  (methanol-water  solution)  and  gas  (carbon  dioxide 
gas,  methanol  vapor  and  water  vapor)  two-phase  mass  transport  in  the  porous  anode  and  cathode  is  formulated  based  on  classical  multiphase  flow 
theory  in  porous  media.  In  the  anode  and  cathode  catalyst  layers,  the  simultaneous  three-phase  (liquid  and  vapor  in  pores  as  well  as  dissolved 
phase  in  the  electrolyte)  water  transport  is  considered  and  the  phase  exchange  of  water  is  modeled  with  finite-rate  interfacial  exchanges  between 
different  phases.  This  model  enables  quantification  of  the  water  flux  corresponding  to  each  of  the  three  water  transport  mechanisms  through  the 
membrane  for  DMFCs,  such  as  diffusion,  electro-osmotic  drag,  and  convection.  Hence,  with  this  model,  the  effects  of  MEA  design  parameters  on 
water  crossover  and  cell  performance  under  various  operating  conditions  can  be  numerically  investigated. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  liquid-feed  direct  methanol  fuel  cell  (DMFC)  has 
recently  received  much  attention  as  it  has  been  regarded  as  an 
attractive  alternative  to  conventional  power  sources  for  portable 
electronic  devices.  However,  the  commercialization  of  DMFC 
technology  is  still  hindered  by  several  technological  problems, 
among  which  water  management  is  one  of  the  key  issues  [1-4]. 
Unlike  in  gas-hydrogen-feed  polymer  electrolyte  fuel  cells 
(PEFCs),  in  DMFCs,  liquid  methanol  solution  is  fed  to  its 
anode.  As  a  result,  liquid  water  is  prone  to  transport  through 
the  membrane  from  the  anode  to  the  cathode,  which  is  usually 
called  water  crossover.  Water  crossover  can  not  only  result  in 
a  water  loss  from  the  anode,  but  also  increase  the  difficulty 
in  avoiding  the  cathode-flooding  problem,  downgrading  the 
DMFC  performance  significantly.  Therefore,  suppressing  water 
crossover  is  beneficial  not  only  for  simplifying  the  DMFC 
system  but  also  for  improving  cell  performance.  To  this  end,  it 
is  essential  to  gain  a  better  understanding  of  the  mechanisms  of 
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water  crossover  in  DMFCs,  which  appears  to  be  significantly 
different  with  that  in  PEFCs  [1-3].  Over  the  past  decade,  the 
problem  of  water  crossover  through  Nafion  membranes  for 
DMFCs  has  been  experimentally  studied  extensively  [1-11]. 
However,  it  is  rather  difficult  to  shed  light  on  each  mechanism 
of  water  crossover  through  the  membrane  through  experimental 
investigations,  as  water  transport  in  DMFCs  is  intrinsically 
coupled  with  a  series  of  physicochemical  processes,  including 
species  transport,  momentum  transport,  and  multiple  electro¬ 
chemical  reactions.  Hence,  mathematical  modeling  plays  an 
important  role,  as  it  can  provide  a  powerful  and  economical  tool 
to  quantify  the  complex  transport  processes  and  thus  elucidate 
the  water  transport  mechanisms  in  DMFCs. 

Recently,  extensive  efforts  have  been  made  to  develop 
two-phase  mass  transport  models  for  DMFCs  [12-17],  which 
are  more  realistic  and  effective  as  the  coexisting  liquid  and 
gas  flow  behavior  has  a  significant  influence  on  the  mass 
transport  processes  [18,19].  Murgia  et  al.  [12]  presented  a 
one-dimensional  model  based  on  phenomenological  transport 
equations  for  DMFCs.  In  order  to  consider  the  two-phase  flow 
interaction  in  the  diffusion  layer  (DL),  they  introduced  a  Gaus¬ 
sian  function  to  approximately  account  for  the  influence  of  the 
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Nomenclature 

water  vapor  activity 

A\g 

interfacial  specific  area  between  liquid  and  gas 
phase  (m2  m-3) 

Ay 

specific  area  (m2  m-3) 

c 

molar  concentration  (molm-3) 

D 

diffusivity  (m2  s-1) 

F 

Faraday  constant  (96,478  C  mol-1) 

h\g 

interfacial  transfer  rate  constant  for  methanol 
(ms-1) 

I 

current  density  (A  m-2) 

h 

parasitic  current  resulting  from  methanol 
crossover  (A  m-2) 

jo 

exchange  current  density  (A  m-2) 

7a 

anode  current  density  (Am-3) 

7  c 

cathode  current  density  (Am-3) 

kc 

condensation  rate  (s-1) 

kQ 

evaporation  rate  (atm  s)-1 

K 

relative  permeability 

K 

permeability  of  porous  material  (m2) 

m 

Source  term  in  mass  conservation  equation 
(kgm-3  s-1) 

M 

molecular  weight  (kg  mol-1) 

n<i 

electro-osmotic  drag  coefficient 

N 

mol  flux  (mol  m-2  s-1),  or  Nation 

Pc 

capillary  pressure  (Pa) 

Pg 

gas  phase  pressure  (Pa) 

P\ 

liquid  phase  pressure  (Pa) 

R 

gas  constant,  J  (mol  K)-1 

R 

source  term  in  species  conservation  equation 
(molm-3  s-1) 

R 

interfacial  species  transfer  rate  (molm-3  s-1) 

^contact 

ohmic  contact  resistance  (£2  m2) 

liquid  saturation 

T 

temperature  (K) 

Vo 

thermodynamic  equilibrium  voltage  (V) 

Pcell 

Cell  voltage  (V) 

X 

coordinate,  m,  or  mole  fraction  in  liquid  solution 

y 

Coordinate,  m,  or  mole  fraction  in  gas  mixture 

Greek  letters 

a 

net  water-transport  coefficient 

«a 

anode  transfer  coefficient  at  anode 

a  c 

cathode  transfer  coefficient  at  cathode 

8 

thickness  of  porous  layer  (m) 

£ 

porosity  of  porous  medium 

y 

reaction  order  of  ORR 

ii 

overpotential  (V) 

K 

ionic  conductivity  of  membrane  (£2-1  m-1) 

Ka(Kd) 

absorption  (desorption)  coefficient  (s-1) 

^adsC^de: 

s)  H20dls  absorption  (H20liq  desorption)  coeffi¬ 
cient  (s-1) 

X 

water  content 

P' 

viscosity  (kgm-1  s-1) 

0  contact  angle  (°) 

p  density  (kgm-3) 

o  interfacial  tension  (N  m-1) 

Superscripts 
eff  effective  value 

in  inlet  condition 

ref  reference  value 

sat  saturated  value 

*  in  equilibrium 

Subscripts 

a  anode 

bf  back  flow 

c  cathode,  or  capillary 

dry  dry  membrane 

dv  from  dissolved  phase  to  vapor  phase 

e  electrolyte,  or  evaporation 

g  gas  phase 

1  liquid  phase 

Id  from  liquid  phase  to  dissolved  phase 
mem  membrane 

M  methanol 

MV  methanol  vapor 

rg  relative  value  for  gas  phase 

rl  relative  value  for  liquid  phase 

vapor  water  vapor 

vl  from  vapor  phase  to  liquid  phase 

W  water 

we  dissolved  water 

WV  water  vapor 


capillary  pressure  on  the  effective  gas  porosity.  Wang  and  Wang 

[13]  modeled  a  DMFC  using  the  multiphase  mixture  model,  in 
which  the  species  in  the  liquid  and  gas  phase  are  assumed  to  be 
at  the  thermodynamic  equilibrium  condition.  Rice  and  Faghri 

[14]  proposed  a  transient,  multiphase  model  for  a  passive  fuel 
feed  DMFC.  Noteworthy  is  that  the  evaporation/condensation 
of  methanol  and  water  was  formulated  in  a  manner  to  capture 
non-equilibrium  effects  between  phases,  which  differed  from 
the  thermodynamic  equilibrium  assumption  employed  else¬ 
where  [12,13].  Most  recently,  Yang  and  Zhao  [16,17]  developed 
a  two-phase  mass  transport  model  for  liquid-feed  DMFCs, 
which  was  formulated  based  on  classical  multiphase  flow  theory 
in  porous  medium  and  with  the  effect  of  non-equilibrium  evap¬ 
oration  and  condensation  of  methanol  and  water  taken  account. 
Most  of  previous  work  about  DMFCs  mainly  investigated  the 
transport  phenomena  in  porous  regions  of  the  anode  and  the 
cathode,  but  relatively  little  efforts  have  been  made  on  modeling 
of  the  simultaneous  water  transport  across  the  membrane. 

Water  transport  across  the  Nation  membrane  has  been 
extensively  modeled  over  the  past  decade  for  hydrogen-feed 
PEFCs  [20-38].  Springer  et  al.  [20]  developed  a  simple,  one¬ 
dimensional  isothermal  model  of  a  complete  PEFC  that  has 
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provided  useful  insight  into  the  cell’s  water  transport  mech¬ 
anism.  They  considered  gas  transport  in  the  electrode  and 
dissolved  water  transport  in  the  membrane,  and  applied  thermo¬ 
dynamics  phase  equilibrium  condition  between  dissolved  water 
in  the  membrane  and  water  vapor  in  the  electrode  at  the  mem¬ 
brane/electrode  interfaces.  This  equilibrium  assumption  has 
been  widely  adopted  by  many  researchers  [21-29]  in  modeling 
the  water  transport  through  the  membrane.  This  model  had  been 
extended  by  Kulikovsky  [30],  who  considered  transport  of  both 
gas  phase  and  dissolved  water  in  electrolyte  phase  in  the  catalyst 
layers  (CLs).  They  also  assumed  that,  in  the  CLs,  the  dissolved 
water  phase  is  in  phase  equilibrium  with  the  water  vapor  in  voids, 
and  thus  water  transport  in  the  CLs  can  be  dealt  with  based 
on  the  “fictitious  water  concentration”  treatment.  This  approach 
can  also  be  found  elsewhere  [25-28].  However,  the  equilibrium 
assumption  between  water  vapor  and  dissolved  water,  albeit 
suitable  as  an  approximation,  is  not  strictly  correct,  since  the 
equilibrium  state  is  not  achieved  instantly  but  over  a  prolonged 
period  of  time  [31].  As  thus,  Siegel  et  al.  [32]  and  Vorobev  et  al. 
[31]  developed  models  for  PEFCs,  in  which  they  did  not  assume 
equilibrium  between  water  vapor  and  dissolved  water  in  the  CLs, 
but  rather  considered  them  as  different  phases  exchanging  at  a 
finite  rate.  Most  recently,  since  under  typical  PEFC  operating 
conditions  both  liquid  water  and  water  vapor  exists  in  voids  of 
the  electrodes,  a  more  realistic  one-dimensional  model  had  been 
developed  by  Shah  et  al.  [33,34],  in  which  three  phases  of  water 
(i.e.,  dissolved  water  in  the  electrolyte,  water  vapor  in  voids  and 
liquid  water  in  voids)  were  considered.  Notably,  the  equilib¬ 
rium  assumption  between  different  phases  was  not  employed, 
and  evaporation  and  condensation  between  liquid  and  vapor,  as 
well  as  adsorption  and  desorption  between  dissolved  water  and 
vapor  and  between  dissolved  water  and  liquid,  were  introduced 
for  the  finite-rate  interfacial  mass  transport.  This  basic  feature 
of  the  model  enables  in  capturing  the  convoluted  water  transport 
phenomena  across  the  whole  MEA  for  PEFCs. 

However,  with  regard  to  DMFCs,  most  of  previous  models 
assumed  that  the  membrane  was  fully  hydrated  and  no  water  dif¬ 
fusion  occurred,  and  little  work  about  the  water  transport  across 
the  membrane  was  reported  [15,39-44].  Meyers  and  Newman 
[39-41]  developed  a  mathematical  model  that  described  the 
transport  of  species  in  a  multi-component  membrane.  Liquid 
phase  in  the  anode  and  gas  phase  in  the  cathode  were  considered, 
and  the  equilibrium  condition  was  assumed  between  the  mem¬ 
brane  and  the  adjacent  liquid  or  vapor  phases.  Similarly,  Schultz 
and  Sundmacher  [42]  developed  a  one-dimensional  rigorous 
process  model  for  a  single-cell  DMFC,  which  also  introduced 
phase  equilibrium  assumption  between  the  dissolved  water  and 
the  water  inside  the  pores  of  both  CLs.  These  two  reports  all 
considered  single  phase  flow  in  both  the  anode  and  the  cathode, 
which  were  just  for  a  first  approximation  and  could  not  reflect 
the  mass  transport  processes  in  real  liquid-gas  two-phase  flow 
occurring  in  the  DMFC.  Recently,  Liu  and  Wang  [15,43]  devel¬ 
oped  a  liquid-gas  two-phase  transport  model  for  DMFCs,  in 
which  the  water-crossover  flux  was  determined  based  on  the  con¬ 
ditions  at  the  CL/membrane  interfaces  and  the  phase  equilibrium 
assumption.  Different  water-crossover  mechanisms,  including 
diffusion,  electro-osmotic  drag  and  convection,  were  analyzed. 


The  objective  of  this  work  is  to  develop  a  one-dimensional, 
isothermal  two-phase  mass  transport  model  for  liquid-feed 
DMFCs.  Water  transport  in  three  phases  (i.e.,  liquid,  vapor,  and 
dissolved  phase)  is  considered  simultaneously,  and  finite-rate 
interfacial  exchanges  between  different  phases  are  introduced 
without  any  phase  equilibrium  assumptions.  With  these  finite- 
rate  interfacial  exchanges  between  different  phases,  dissolved 
water  transport  across  the  membrane  can  be  intrinsically  cou¬ 
pled  with  the  liquid-gas  two-phase  flow  in  the  porous  regions 
of  both  the  anode  and  the  cathode.  Thus,  this  model  enables 
us  to  investigate  different  water-crossover  mechanisms  sepa¬ 
rately,  and  study  the  effect  of  cell  operating  conditions  on  water 
crossover.  Furthermore,  since  the  cathode  porous  region,  espe¬ 
cially  the  cathode  micro-porous  layer  (MPL),  plays  an  important 
role  in  both  water  crossover  through  the  membrane  and  water 
ejection  from  the  cathode  [1,2,45-47],  this  model  also  includes 
MPLs  between  the  DL  and  CL  of  both  the  anode  and  the  cath¬ 
ode,  which  enables  us  to  study  the  effect  of  membrane  electrode 
assembly  (MEA)  structure  on  water  transport  in  the  DMFC.  The 
results  presented  in  this  work  provide  a  useful  insight  into  the 
water-crossover  mechanisms,  which  is  helpful  for  the  optimal 
design  of  the  MEA  to  reduce  water  crossover  and  to  achieve 
better  cell  performance. 

2.  Formulation 

We  consider  a  MEA,  as  illustrated  in  Fig.  1 ,  which  consists  of 
seven  layers  from  the  interface  I  between  the  anode  channel  and 
anode  diffusion  layer  (ADL)  to  the  interface  VIII  between  the 
cathode  diffusion  layer  (CDL)  and  cathode  gas  channel.  Note 
that  micro-porous  layers  (MPLs)  are  also  included  between  the 
catalyst  layers  (CLs)  and  coarse  DLs  at  both  the  anode  and 
the  cathode.  At  the  anode,  liquid  methanol  solution  is  fed  from 
the  channel  through  the  DL  to  the  CL,  while  the  produced  gas 
CO2  is  expelled  from  the  DL  to  the  channel.  Thus,  a  liquid-gas 
two-phase  counter  flow  is  involved  in  the  anode  porous  region, 


I  *  II  III  V  VI  VII  VIII 


Fig.  1.  Schematic  of  the  model  domain. 
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Fig.  2.  Schematic  representation  of  the  interaction  between  water  in  the  three 
phases  in  the  catalyst  layers. 


including  the  ADL,  the  AMPL  and  the  ACL.  Similarly,  in  the 
gas-feed  cathode,  liquid  water  usually  resides  in  the  porous 
region  under  typical  operating  conditions,  and  thus  a  liquid-gas 
two-phase  counter  flow  also  occurs.  In  the  following,  we  present 
a  one-dimensional  mass  transport  model  for  water  through  the 
MEA. 


2.7.  Phase  changes  of  water  in  the  CLs 

In  the  CLs,  both  liquid  and  water  vapor  can  dissolve  in  the 
electrolyte  (typically  Nafion)  with  the  effect  of  adsorption  to 
form  dissolved  water.  Thus,  as  schematically  illustrated  in  Fig.  2, 
water  exists  in  three  phases  in  the  CLs:  liquid,  vapor,  and  the 
dissolved  phase.  To  model  the  water  transport  in  the  CLs,  the 
phase-change  rates  between  different  phases  are  needed.  As 
shown  in  Fig.  2,  both  the  water  produced  at  the  cathode  and 
that  consumed  at  the  anode  is  in  the  liquid  phase.  For  the  phase 
change  between  liquid  and  vapor,  the  rate  of  condensation  and 
evaporation  can  be  modeled  using  the  finite-rate  approach  [17]: 

*vl  =  hY\(x- vapor Pg  ~  P vapor)  (1) 


In  a  similar  fashion,  the  dissolved  water-vapor  transfer  rate 
through  desorption  and  adsorption  can  be  expressed  as: 


7?dv  —  ^dvC^we  CWe  v) 


(3) 


where  Cwe  is  the  dissolved  water  concentration  in  the  electrolyte 
and  C^e  v  is  the  equilibrium  dissolved  concentration  when  the 
polymer  is  in  equilibrium  with  water- vapor- saturated  gas.  The 
driving  force  (Cwe  —  C^e  v)  is  positive  for  desorption  but  nega¬ 
tive  for  adsorption.  The  dissolved  water  concentration  Cwe  can 
be  transformed  to  the  water  content  A  in  the  electrolyte  and  the 
relationship  between  Cwe  and  A  is  given  by 

EW 

A  =  Cwe  (4) 

Pdry 

In  Eq.  (3),  C^e  v  can  be  derived  from  the  equilibrium  water 
content  A*  e>v  at  353  K  [48]  to  give 

^we,v  @353  k  =  0-3  +  10.8aw  —  16  a^  +  14.1  a^  (5) 

and  at  303  K  [20] 

^we,v@303K  =  0-043  +  17.81aw  —  39.85 +  36. 0a^  (6) 

where  aw  is  the  water  vapor  activity  (aw  =  vVapor pg/ Pvapor)-  The 
equilibrium  water  content  A^e  at  any  temperature  can  thus  be 
approximated  by  [49] 


^we,v  ^we.v  @303K  ~h 


^we,v  @353  K  ^we,v  @303  K 

50 


(T  -  303) 


(7) 


The  coefficients  of  adsorption  and  desorption  h& v  can  be  given 
by  [33,34] 


/?dv 


1 

-/cd(l  -  s) A 


+ 


|C* 


where  /4apor  *s  the  saturation  pressure  of  water  vapor,  pg  is  the 
gas  pressure,  and  vVapor  is  the  molar  fraction  of  water  vapor  in  the 
gas  phase.  Eq.  (1)  indicates  that  the  phase  change  between  liquid 
and  vapor  is  driven  by  the  deviation  of  the  local  thermodynamic 
state  from  equilibrium.  Note  that  the  driving  force  (xVaporPg  — 
7* vapor)  is  positive  for  condensation  but  negative  for  evaporation. 
The  mass-transfer  coefficient  hy\  can  be  given  by  [17,33]: 

1  l-^vaporPg  —  Pvaporl 
*  vapor  Pg  —  Pvapor 

_l_  ke£Sp\  (  kvaporPg  “  Pvaporl 

2Mh20  y  3C  vapor  Pg  —  P  vapor 


Avi  — 


kc£(  1  —  s)x 


vapor 


2  RT 


+  ^a(l  ~s)k 


1 Cwe 

CWe 


(8) 


where  k&  and  /ca  are  the  desorption  and  adsorption  rate  constants. 

It  is  known  that  when  the  electrolyte  Nafion  is  submerged 
in  liquid  water  its  equilibrium  water  content  appears  to  jump 
discontinuously  to  a  higher  value  of  A^e  x  —  22  [34,39-41]. 
To  capture  this  anomaly  we  also  introduce  the  liquid-dissolved 
transfer  rate  through  desorption  and  adsorption,  which  is  driven 
by  the  deviation  of  the  water  concentration  from  the  equilibrium 
one  and  is  given  by 


^ld  —  ^ld(Cwe,i  Cwe) 


(9) 


where  kc  and  kQ  are  the  condensation  and  evaporation  rate  con¬ 
stants,  £  the  porosity  of  the  porous  region,  s  the  liquid  water 
saturation,  T  the  cell  temperature,  p\  the  density  of  liquid  and 
Mh2o  the  molecular  weight  of  water.  In  the  above  expression,  a 
sufficiently  high  values  of  kc  and  Ae  are  chosen  so  that  the  vapor 
is  essentially  in  equilibrium  with  the  liquid. 


where  C^e  j  is  the  equilibrium  dissolved  concentration  when  the 
electrolyte  is  in  equilibrium  with  liquid  water.  The  driving  force 
(C*e  j  —  Cwe)  is  positive  for  desorption  of  water  from  the  elec¬ 
trolyte  to  form  liquid  but  negative  for  adsorption  of  liquid  water 
to  the  electrolyte  in  the  CLs.  The  coefficients  of  adsorption  and 
desorption  h\&  is  assumed  to  be  dependent  of  liquid  saturation 
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and  can  be  given  by  [33,34] 


Aid 


1 

—  2^adsS 


|Cwe 

Cwe 


~  ^wejlA 

-  qje.1  J 


1 

+  -/CdesS 


I  Cwe 

^we 


(10) 


where  k des  and  /ca ds  are  the  desorption  and  adsorption  rate  con¬ 
stants  between  liquid  and  dissolved  water.  The  values  of  k des 
and  /cads  are  chosen  large  enough  to  avoid  significant  overshoot 
of  cwe  when  it  exceeds  C^e  These  parameters  and  their  values 
are  given  in  Table  1 . 


2.2.  Governing  equations  for  mass  transport 


We  now  present  the  governing  equations  for  the  one¬ 
dimensional,  steady-state,  isothermal  two-phase  mass  transport 
in  the  porous  regions  of  the  DMFC  anode  and  cathode,  as  well  as 
for  the  dissolved  water  transport  in  the  electrolyte  in  both  the  CLs 
and  the  membrane.  The  details  of  the  two-phase  mass  transport 
model  in  the  porous  regions  can  be  found  elsewhere  [16,17]. 
Referring  to  Fig.  1,  in  the  anode  and  cathode  porous  regions 
(from  Interface  I  to  IV  and  from  V  to  VIII),  there  are  a  total  of 
nine  variables  that  are  used  to  describe  mass  transport  of  each 
species,  including  the  concentration  of  liquid  methanol  (Cm,i), 
the  concentration  of  methanol  vapor  (CMV,g),  gas  pressure  in 
the  anode  (/?g?a),  liquid  pressure  in  the  anode(/?i?a),  the  concen¬ 
tration  of  water  vapor  in  the  anode  (Cwv,g,a),  the  concentration 
of  gas  oxygen  in  the  cathode  (Co2,g),  liquid  phase  pressure  in 
the  cathode  (pi?c),  gas  phase  pressure  in  the  cathode  ( pg?c )  and 
the  concentration  of  water  vapor  in  the  cathode  (Cwv,g,c)«  In  the 
electrolyte  region  (from  III  to  VI),  variable  of  dissolved  water 
concentration  in  the  electrolyte  (Cwe)  are  involved.  The  govern¬ 
ing  equations  that  describe  the  mass  conservation  of  each  species 
at  different  phases  are  given  below: 


Anode  porous  region  (I-IV): 


Pi, a  : 

v-( 

Kkv\  „  \ 

Ml/Pl  PU) 

=  mu 

(id 

Rg,a  : 

v( 

<  Mg/ Pg 

)  =  '«g,a 

(12) 

Cm,i  • 

V  . 

)  Cm,i  -  VCm,i 

=  ^M,l,a 

(13) 

Kkvg 

Pg 


CMV,g  —  D 


eff 

MV,g 


VCMv,g 


(14) 


Cwv,g  :  V  • 

=  %V,g,a 


6wV,g,a 


—  ^WY,g^^WV,g,a 


(15) 


Cathode  porous  region  (V-VIII): 


Pg,c 


Pl,c 


V  • 


V  • 


Kh 


rg 


P  g/ Pg 
KkvX 
P\/p\ 


Vpg,c 

Vpi,c 


m 


g.c 


mu 


Co2,g  •  V  • 

=  ^02, g 


Kh 


rg. 


Pg 


V/?g,c  1  Co2,g  ^02,gV^02,] 


(16) 

(17) 


(IB) 


Cwv,g,c  •  V 

=  %V,g,c 


Kh 


Pg 


rgV/7g,c  )  Cwv,g,c  -  7)wv,gV(7wv,g,c 


(19) 


Here,  it  should  be  noted  that  the  capillary  pressure  for  the 
porous  region  is  given  by  [16,17]: 

pc  =  pg  —  px  =  a  cos  0(s/K)°'5  J(s)  (20) 

where  J(s)  represents  the  widely-used  Leverette  function: 


J  1.417(1  -  s)  -  2.120(1  -  s)2  +  1.263(1  -  s)3  0  <  0  <  90° 

J(S)  ~  \  1.4175  -  2.12052  +  1.2635s  90  °  <0<  180° 

(21) 

Thus,  the  liquid  saturation  for  the  porous  regions  of  both 
the  anode  and  the  cathode  can  be  determined  from  Eqs. 
(20)  and  (21)  after  obtaining  the  liquid  pressure  and  gas 
pressure. 


2.3.  Electrolyte  region  (III-VI) 


Unlike  in  the  porous  region,  in  the  electrolyte  region  only 
dissolved  water  and  methanol  need  to  be  considered,  as  the 
membrane  is  usually  regarded  as  a  gas  insulator.  Transport  of 
dissolved  water  through  the  electrolyte  generally  depends  on 
molecular  diffusion,  electro-osmotic  drag  and  back  convection. 
Accordingly,  the  flux  of  water  crossover  (Vw)  can  be  given  by 
[1,2]: 

NW  =  -e‘-5Dwea)VCwe  +  nd,H2o7  -  £e^emAVpi  (22) 

F  p\Mu2o 

Thus,  the  governing  equation  for  the  dissolved  water  concen¬ 
tration  (Cwe)  is: 


V  •  Vw  =  V  •  (-£g-5Z)We(A.)VCWe)  +  V  •  (  Hd,H20  — 


Rid  —  Rdv  —  Rbf(ACL) 
0(PEM) 

Rid  —  Rdv  +  Rbf(CCL) 


(23) 


Here,  we  assume  that  the  water-crossover  flux  due  to  back  con¬ 
vection  (—  V p\)  is  uniformly  distributed  across  the  CLs, 
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Table  1 

Physicochemical  properties 


Parameters 

Symbols 

Value 

Unit 

Ref. 

Porosity,  permeability 

ADL 

£adl  Wadi 

0.75,  1.0  x  10-12 

m2 

- 

AMPL 

£ampl  Warnpl 

0.3, 2.5  X  10“13 

m2 

[16] 

ACL 

^acl  Wacl 

0.3,  1.5  x  10“ 14 

-,  m2 

[16] 

MEM 

£mem  Wrnem 

0.3,  2.0  x  10“18 

-,  m2 

[16] 

CCL 

£ccl  ,*ccl 

0.3,  1.5  x  10“ 14 

-,  m2 

[16] 

CMPL 

£cmpl  Wcmpl 

0.3, 2.5  x  10-13 

-,  m2 

[16] 

CDL 

£cdl  Wcdl 

0.75,  1.0  x  10“12 

-,  m2 

- 

Nation  volume  fraction 

ACL 

£e,acl 

0.3 

- 

- 

CCL 

£e,ccl 

0.3 

- 

- 

Diffusivities 

MeOH  in  water 

Dm,i 

1.58  X  10 -9e0.02623(r-298) 

m2  s  1 

[16] 

MeOH  in  Nation 

Dm,n 

4.9  x  1 0“  10e[2436O/333-  1/;T)] 

m2s-1 

[16] 

Methanol  vapor 

DM,g 

-6.954  X  10“6 +  4.5986  x  10_8r+ 9.4979  x  10“ur2 

m2s-1 

[17] 

O2  in  gas 

°02,g 

1.775  x  10-5  (575+5)  * 

m2  s_1 

[16] 

Water  vapor 

Dwv,g 

2.56  x  10— 5  (  307715  )  2  334 

m2  s-1 

[16] 

Dissolved  water  in  Nation 

12  wc 

4.17  x  10_8A(161e_)l  +  l)e_2436/r 

m2  s_1 

[21] 

Viscosity  of  gas  phase 

Mg 

2.03  x  10“5 

kgm-1  s_1 

[17] 

Viscosity  of  liquid  phase 

Mi 

4.05  x  10"4 

kgm-1  s_1 

[16] 

Electro-osmotic  drag  coefficients  of 

^d,H20 

2.5, 

22  A 

- 

[24] 

water  and  methanol 

^d,M 

«d,H20*M 

- 

- 

Evaporation  rate  constant  for  water 

1.0  X  103 

(atms)-1 

- 

Condensation  rate  constant  for  water 

kc 

1.0  X  106 

s-1 

- 

Interfacial  transfer  rate  constant  for 

h\ g 

0.05 

ms-1 

- 

methanol 

Specific  interfacial  area  between 

Aig 

105 

m-1 

[17] 

liquid  and  gas 

Water  vapor  adsorption  (desorption) 

KaOd) 

8.0  (0.25) 

s  1 

- 

coefficient 

Water  liquid  adsorption  (desorption) 

^ads(^des) 

4.0  (100.0) 

s  1 

[34] 

coefficient 

Proton  conductivity  in  membrane 

K 

7  3e[1268(l/298-l/r)] 

[16] 

Henry  law  constant  for  methanol 

kn,M 

0.096e°-°4511(T_273) 

atm 

[17] 

The  saturation  pressure  of  water 

l°gl0PwY 

-2.1794  +  0.02953(7-273)  -9.1837  x  10-5(T-273)2  +  1.4454  x  10-7(T-273)3 

atm 

[17] 

vapor 

The  saturation  pressure  of  methanol 

Pmy 

kH*M,l 

atm 

[17] 

vapor 

Thermodynamic  voltage 

Vo 

1.21 

V 

[16] 

Transfer  coefficient  of  anode 

aa 

0.5 

- 

[16] 

Transfer  coefficient  of  cathode 

occ 

1.0 

- 

[16] 

Anode  exchange  current  density 

A  /ref 

1.0  x  105 

Am-3 

[16] 

Cathode  exchange  current  density 

A  ;ref 
/iv>c7o,02 

6.97  x  102 

Am-3 

[16] 

Anode  reference  concentration 

z^-ref 

100 

mol  m-3 

[16] 

Cathode  reference  concentration 

/^ref 

C02 

36.5 

mol  m-3 

[16] 

Surface  tension 

cr 

0.0644 

N  m_1 

[45] 

Equivalent  weight  of  ionomer 

EW 

1.1 

kg  mol-1 

[21] 

Dry  membrane  density 

Pdry 

1980 

kgm-3 

[21] 

and  the  term 


Rbf  = 


M^h2o  ^mem  <$cl 


(24) 


is  to  account  for  the  effect  of  back  convection,  where  A/?ijC_a 
represents  the  liquid  pressure  difference  between  the  cathode 
and  the  anode.  The  flux  of  methanol  crossover  (A/m,i),  which 
depends  on  the  effect  of  molecular  diffusion,  electro-osmotic 


drag  and  convection,  can  be  given  by: 


/ 

Nm,\=  —  T>m,nVCm,1  +  —  - 


(  ^mem  A/?i?c_a\ 

- I -  Cm’! 

\  M7  Omem  / 


(25) 


It  should  be  mentioned  here  that  the  model  developed  in  this 
work  can  also  be  used  to  study  the  effect  of  methanol  crossover 
on  cell  performance. 

To  make  the  above  governing  equations  closed,  some  con¬ 
stitutive  correlations  and  definitions  are  needed.  These  include 
relative  permeability  for  both  gas  and  liquid  phases,  effective  dif- 


Table  2 

Constitutive  relations 
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Parameters 


Relative  permeabilities 


Effective  diffusion  coefficients  of  species  [16,17] 


General  generation  rate  of  mass  in  liquid  phase 


General  generation  rate  of  mass  in  gas  phase 


Mole  generation  rate  of  species 


Expressions 


kr\=s3 
krg  =  (l-s)3 

Dfg  =  A,g£L5(l  -  s)L5  i  :  02,  WV,  MV 

DM,\sh5sh5 
(e  +  £e) 

[sADmjs'V-5)  +  Se/(CM,N£eL5)] 
D  m,n 


m,  a  = 


/^H20  Al  —  %^MV,g 

Mh2o(Ai  -  Ad  +  Af  -  ja/6F)  -  Mm(R MV,g  +  ja/6F  +  Ip/6F8ac\ ) 


= 


f 

|  ^H2' 


O  Al 

o(Al  —  Ad  —  Af  Hb  ((7c  —  /p/<5cci)/2F  +  Ip/3F8CC\)) 


_  I  -Mr2oAi  +  /?MV,g 

m2  a  Mco2  ^76^  +  Mh2o(4  -  Ai)  +  Atv,g 


-Mh2oAi 

-Mq2  7c /4F  +  Mco2  7p/6F(5cci  +  Mh2o(Av  -  Al) 


f  0  •  f 

S  7  c  >%V,g,c  =  < 

l  4E  l 


— Ai 
Av  —  Ai 


A2,g  =  ^  7c  ,Rwv,g,c  = 

4F 

Av,g 

7a  •  /p  »  %V,g,a  = 

_6F  "  RmV,s  “  6Riaci 

D  4  7/1  \(%y_^MV) 

^MV,g  —  AighigS\  1  5') 


— Ai 
Av  —  Al 


Liquid 

Gas 

ADL/AMPL  ACL  MEM 


ADL/AMPL  ACL 

CDL/CMPL  CCL 

ADL/AMPL  ACL 

CDL/CMPL  CCL 

CDL/CMPL  CCL 

ADL/AMPL  ACL 

ADL/AMPL/ ACL 


fusion  coefficients  for  each  species,  the  interfacial  transfer  rates 
of  methanol  between  liquid  and  gas  phases  and  all  the  other 
source  terms.  All  these  correlations  and  associated  nomencla¬ 
tures  are  listed  in  Table  2. 

2.4.  Boundary  and  interfacial  conditions  for  mass  transport 

As  indicated  in  Fig.  1,  the  computational  domain  is  enclosed 
by  eight  boundaries/interfaces.  The  conditions  at  each  bound¬ 
ary/interface  are  described  below. 

Boundary  I:  This  boundary  represents  the  inlet  of  reactant  sup¬ 
ply  at  the  anode,  at  which  the  concentration  of  liquid  methanol, 
liquid-phase  pressure,  liquid  saturation  and  the  concentrations 
of  methanol  vapor  and  water  vapor  are  all  specified  to  be  inlet 
conditions: 


Interface  IV:  This  interface  is  the  interface  between  the  ACL 
and  the  membrane  for  the  anode,  which  is  impermeable  for 
both  the  gas  and  the  liquid  (except  methanol  in  the  liquid). 
Thus,  all  the  fluxes  except  the  dissolved  water  and  methanol  in 
the  liquid  in  the  x  direction  are  zero  at  this  interface: 


dpi, a. 

_  q  dPg,a 

_  q  3CwV,g,a 

Q 9CMV,g 

dx 

dx 

dx 

dx 

(28) 

Interface  V:  This  interface  is  the  interface  between  the  CCL  and 
the  membrane  for  the  cathode.  It  is  assumed  that  the  methanol 
transported  from  the  anode  will  be  depleted  immediately  due 
to  the  very  fast  methanol  reaction  at  the  cathode,  and  thus  the 
methanol  concentration  at  this  interface  is  zero.  Similar  to  inter¬ 
face  IV,  all  the  gas  and  liquid  fluxes  in  the  x  direction  are  zero 
at  this  interface: 


Cm,\  —  ^M’  Cmv  —  cjjjy,  Cwv  —  ^wv,a ,  Pi  —  Pi,  s  —  1 

(26) 

Interface  III  and  VI:  These  two  boundaries  are  the  interfaces 
between  MPLs  and  CLs  for  the  anode  and  cathode,  which 
indicate  impermeable  walls  for  the  dissolved  water  since  there 
exists  no  electrolyte  in  the  MPLs.  Accordingly,  the  fluxes  of  the 
dissolved  water  in  the  x  direction  are  zero  at  these  interfaces: 


I  =0,^ 

0 dPz* 

3Cwv,g,c 

1  dx 

+  ’  dx 

+  ’  dx 

Cm,\\_  =  1 


=  0, 


3Cq2 


dx 


=  0 


(29) 


Boundary  VIII:  This  boundary  represents  the  inlet  of  oxy¬ 
gen  supply  at  the  cathode,  at  which  the  following  boundary 
conditions  are  specified: 


(27)  Cq2  =  C$2,  Cwv  =  Cwv.c’  Pg  =  Pg'  s  =  0 


Nw\+  =  0(111),  Nw\-=  0(VI) 


(30) 
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Table  3 


Cell  geometric  dimensions  and  operating  parameters 


Parameters 

Symbols 

Value 

Unit 

Anode  diffusion  layer  thickness 

<$adl 

2.6  x  10“4 

m 

Anode  micro-porous  layer  thickness 

^ampl 

0.3  x  10"4 

m 

Anode  catalyst  layer  thickness 

^acl 

0.2  x  10“4 

m 

Membrane  thickness  (Nation  115) 

^mem 

0.5  x  10“4 

m 

Cathode  diffusion  layer  thickness 

<$cdl 

2.6  x  10"4 

m 

Cathode  micro-porous  layer  thickness 

^cmpl 

0.3  x  10“4 

m 

Cathode  catalyst  layer  thickness 

<$ccl 

0.2  x  10"4 

m 

Operation  temperature 

T 

333.15 

K 

Anode  inlet  pressure 

pT 

1.013  x  105 

Pa 

Cathode  inlet  pressure 

pt 

1.013  x  105 

Pa 

Inlet  methanol  concentration  at  anode 

s~<va 

'“'M 

1000 

mol  m" 

-3 

Inlet  methanol  vapor  concentration  at  anode 

/">in 

lMV 

/^sat 

'“'MV 

mol  in' 

-3 

Inlet  oxygen  concentration  at  cathode 

con2 

7.68 

mol  in' 

-3 

Inlet  water  vapor  concentration  at  anode 

/^in 

LWV,a 

P\%oJRT 

mol  m" 

-3 

Inlet  water  vapor  concentration  at  cathode 

/^>in 

LWV,c 

Pl%oJRT 

mol  m" 

-3 

Inlet  liquid  saturation  at  anode 

5in 

1 

- 

Inlet  liquid  saturation  at  cathode 

Sin 

0 

- 

Besides,  the  conditions  at  the  interfaces  II,  III,  VI  and  VII  for 
the  gas  and  liquid  transport  in  the  porous  region  and  at  interfaces 
IV  and  V  for  the  dissolved  water  transport  in  the  electrolyte  are 
given  based  on  the  principle  that  the  continuity  and  mass/species 
flux  balance  are  required  at  each  interface  to  satisfy  the  general 
mass  and  species  conservation  of  the  entire  cell.  It  should  be 
pointed  out  here  that  since  the  anode  gas  and  liquid  pressure,  as 
well  as  the  anode  capillary  pressure,  is  continuous  at  interfaces 
II  and  III,  so  is  the  cathode  gas  and  liquid  pressure  at  interfaces 
VI  and  VII,  the  difference  in  the  characteristics  of  DL,  MPL  and 
CL  causes  discontinuities  in  the  liquid  saturation  at  these  four 
interfaces. 


by 

/=  [  Jadv  (33) 

J  ACL 

To  account  for  methanol  crossover,  the  ‘parasitic’  current 
density  is  used  to  express  the  rate  of  methanol  crossover: 

/p  =  6T7Vm,i  (34) 

where  the  molar  flux  of  methanol  crossover,  Nm,\,  is  given  by 
Eq.  (25). 

On  the  cathode,  it  is  assumed  that  both  the  cell  current  and 
the  ‘parasitic’  current  are  entirely  consumed  by  the  ORR,  i.e.: 


2.5.  Electrochemical  kinetics 

On  the  DMFC  anode,  the  Tafel-like  expression  is  used  to 
model  the  kinetics  of  methanol  oxidation  reaction  (MOR): 

.  .ref  (  CM  V  faaF  \ 

7a  =  ^v,a7o,MeOH  I  I  exP  J  (31> 

where  the  reaction  order  y  is  related  to  the  methanol  con¬ 
centration  and  assumed  to  be  the  zero-order  when  methanol 
concentration  is  higher  than  a  reference  value.  Otherwise,  the 
first-order  reaction  is  specified. 

With  respect  to  the  oxygen  reduction  reaction  (ORR)  on  the 
cathode,  the  Tafel-like  expression  gives  [13,15]: 

7c  =  (1  -  ^)'4v.c70eO2  j  eXp  (_^7?7c)  (32) 

2.6.  Current  balance,  cell  voltage  and  water  balance 

The  protons  and  electrons  produced  by  the  MOR  on  the  anode 
transfer  to  the  cathode  through  the  membrane  and  the  external 
circuit,  respectively.  The  cell  current  density  can  be  calculated 


7  +  /p=/  jcdx  (35) 

J  CCL 

In  summary,  for  a  given  anode  overpotential,  ??a,  the  cell 
current  density,  7,  and  the  ‘parasitic’  current  density,  7p,  can 
be  determined  from  Eqs.  (33)  and  (34),  respectively.  Then,  the 
cathode  mixed  overpotential,  ??c,  with  the  effect  of  methanol 
crossover,  can  be  obtained  from  Eqs.  (32)  and  (35).  Finally,  the 
cell  voltage  can  be  determined  from: 

Kc„«,  +  ^)  (36) 

where  Vo,  ^Contact  and  k  denote  the  thermodynamic  equilib¬ 
rium  voltage  of  a  DMFC,  the  contact  resistance  and  the  proton 
conductivity  of  the  membrane,  respectively. 

Finally,  the  flux  of  water  crossover  with  the  effect  of  diffusion 
(Nw,diffX  electro-osmotic  drag  (A^w,eo)  is  determined  using  the 
values  at  the  interface  IV: 


Vcell  =  Vo  -  rfo  +  ric-  I 


Vw,diff —  £e  7)We(^)VCwe|  +  (37) 

7 

Nw,eo  =  +nd,n2o—  (38) 

t  + 

And  the  total  flux  of  water  crossover  is  also  represented  by 
the  so-called  net  water- transport  coefficient,  a ,  which  is  given 
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by: 


The  above-described  governing  equations  for  the  cell  geo¬ 
metric  dimensions  and  operating  parameters  listed  in  Table  3 
subjected  to  electrochemical  properties  listed  in  Table  1, 
were  solved  numerically  using  a  self-written  code,  which  was 
developed  based  on  the  SIMPLE  algorithm  with  the  finite- 
volume-method  [16,17]. 

3.  Results  and  discussion 

3.1.  Distributions  and  various  species  at  different  current 
densities 

This  section  presents  the  distributions  of  several  species, 
including  concentration  of  methanol  solution,  cathode  oxygen 
concentration,  anode  gas  pressure  and  liquid  saturation,  cathode 
liquid  pressure  and  liquid  saturation,  and  the  water  content  in 
the  electrolyte  region,  at  different  current  densities. 

Fig.  3a  shows  the  distribution  of  methanol  concentration  in 
the  DMFC  anode  and  the  membrane  at  different  current  densi¬ 
ties.  It  can  be  seen  that  at  a  given  current  density  the  methanol 
concentration  decreases  in  each  layer  and  reaches  zero  at  the 
PEM/CCL  interface.  The  different  slopes  of  methanol  concen¬ 
tration  across  different  porous  layers  (ADL,  AMPL,  ACL,  and 
PEM)  indicate  the  different  effective  diffusivities  of  methanol 
(D^[)  in  different  layers.  Fig.  3a  also  shows  that  with  the 
increase  in  current  density,  methanol  concentration  decreases  in 
all  the  layers  as  more  methanol  is  consumed  in  the  ACL.  When 
the  current  density  reaches  about  484.0  mA  cm-2  (^a  =  0.6  V), 
methanol  concentration  approaches  almost  zero  in  the  ACL, 
meaning  that  methanol  is  almost  completely  consumed  in  the 
ACL  and  correspondingly,  484.0  mA  cm-2  is  just  the  methanol- 
transport-controlled  limiting  current  density. 

The  distribution  of  gas  oxygen  concentration  in  the  cathode 
porous  region  at  different  current  densities  is  shown  in  Fig.  3b. 


Distance  across  the  ADL/AMPL/ACL/PEM  (gm) 


Oxygen  is  transported  through  the  gas  pores  from  the  cathode 
channel  to  the  CCL,  in  which  the  oxygen  concentration  becomes 
lower  due  to  the  consumption  of  oxygen  by  the  ORR.  With 
the  increase  in  current  density,  the  oxygen  concentration  is  also 
reduced  due  to  the  increased  oxygen  consumption.  However,  the 
decrease  in  oxygen  concentration  from  the  channel  to  the  CL  is 
relatively  small  even  at  the  largest  current  density,  indicating 
that  the  mass  transport  resistance  of  oxygen  is  rather  low.  This 
is  due  to  faster  transport  of  oxygen  in  the  gas  phase  and  lower 
liquid  water  saturation  in  the  cathode  porous  region,  as  will  be 
discussed  later. 

The  distribution  of  gas  pressure  in  the  anode  at  different 
current  densities  is  shown  in  Fig.  4a.  Due  to  the  hydrophobic 
characteristic  of  the  anode  porous  region,  the  gas  pressure  is 
smaller  than  the  anode  liquid  pressure,  which  is  nearly  a  con¬ 
stant  of  101.3  kPa.  It  can  be  seen  that  the  gas  pressure  can  vary 
greatly  from  the  anode  channel  to  the  ACL.  Apparently,  the  gas 
pressure  in  the  region  close  to  the  PEM  is  higher  than  that  in 
the  region  close  to  the  anode  flow  channel,  as  gas  CO2,  gen¬ 
erated  by  MOR  in  the  ACL,  is  transferred  from  the  ACL  to 
the  flow  channel.  Across  the  interfaces  II  and  III  the  gas  pres¬ 
sure  is  continuous  while  the  gradients  are  different  due  to  the 
different  properties  ( K  and  krg)  at  different  layers,  which  can 
be  understood  by  referring  to  Eq.  (12).  With  the  increase  in 
current  density,  the  gas  pressure  increases  due  to  the  increased 
generation  of  gas  CO2.  It  should  be  pointed  out  here  that  in 
the  ACL  the  gas  pressure  does  not  always  decrease  from  the 
interface  IV  to  III,  and  this  anomalous  distribution  is  caused 
by  the  complicated  mass  exchanges  between  different  phases 
(i.e.,  gas,  liquid,  dissolved  phase)  in  the  CL.  By  referring  to  Eq. 
(12)  and  the  corresponding  source  terms  for  the  gas  transport, 
there  exists  not  only  the  CO2  generation  and  water  vapor  gaining 
from  the  desorption  of  dissolved  water  into  the  gas  phase,  but 
also  the  water  vapor  loss  by  condensing  into  the  liquid  phase  in 
the  CL.  The  fluxes  of  water  exchange  depend  highly  on  the  local 
water  status  (e.g.,  concentrations  of  water  vapor  and  dissolved 
water),  which  in  turn  affects  the  anomalous  distribution  of  gas 
pressure. 


Fig.  3.  Distributions  of  (a)  methanol  concentration  in  liquid  phase  through  the  anode  electrode  and  membrane,  and  (b)  gas  oxygen  concentration  through  the  cathode 
at  different  current  densities. 
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(a) 


Fig.  4.  Distributions  of  (a)  gas  pressure  in  the  anode,  and  (b)  liquid  saturation  in  the  anode  at  different  current  densities. 


The  corresponding  distribution  of  anode  liquid  saturation  at 
different  current  densities  is  shown  in  Fig.  4b.  It  is  seen  that 
anode  liquid  saturation  is  always  smaller  than  1.0  in  the  porous 
region  due  to  the  generation  of  CO2 .  It  deceases  from  the  channel 
across  the  ADL  as  the  liquid  is  transported  from  the  channel  to 
the  CL.  It  is  clear  that  there  is  a  discontinuity  of  anode  liquid 
saturation  at  the  interface  II,  which  is  caused  by  the  continuity  of 
capillary  pressure  across  the  interface.  At  the  interface  II,  from 
the  definition  of  capillary  pressure  (Eq.  (20)),  it  gives: 

cos  Oadl  ^  ^DL  ^  J (aadl,ii) 

/  \  0.5 

=  COS  0AMPL  (  „AMPL  )  J(s AMPL,Il)  (40) 

V  A  AMPL  j 

Thus,  the  large  differences  in  permeability  and  porosity 
between  the  ADL  and  AMPL  cause  a  discontinuity  in  the  liquid 
saturation  across  the  interface  of  these  two  layers.  Since  the  per¬ 
meability  of  the  MPL  is  extremely  small,  the  liquid  saturation 
in  this  finer  porous  region  is  rather  low  (~0. 6-0.7).  Similarly, 
due  to  the  large  difference  in  permeability  and  contact  angle 


between  the  AMPL  and  ACL,  a  discontinuity  in  the  liquid  sat¬ 
uration  occurs  across  the  interface  III,  and  the  less  hydrophobic 
property  ( 0 )  of  the  ACL  mainly  contributes  to  the  relatively  high 
liquid  saturation  (~0.85-0.95).  With  the  increase  in  current  den¬ 
sity,  the  anode  liquid  saturation  is  reduced  due  to  the  increased 
generation  of  gas  CO2.  Generally,  the  gas  phase  in  the  anode 
porous  region  has  two  effects:  on  one  hand,  it  limits  the  access 
of  methanol  from  the  channel  to  the  catalyst  region;  on  the  other 
hand,  the  gas  occupancy  in  the  catalyst  layer  reduces  the  water 
content  in  the  electrolyte  and  thus  influences  the  water  crossover 
from  the  anode  to  the  cathode,  as  will  be  discussed  later. 

The  gas  pressure  in  the  cathode  is  also  nearly  a  constant  of 
101.3  kPa,  while  the  liquid  pressure  can  vary  greatly  from  the 
cathode  channel  to  the  CCL,  which  is  shown  in  Fig.  5a.  Due 
to  the  hydrophobic  characteristic  of  the  cathode  porous  region, 
the  liquid  pressure  is  larger  than  the  cathode  gas  pressure.  It  is 
seen  that  the  liquid  pressure  increases  from  the  cathode  channel 
to  the  CCL  due  to  the  fact  that  liquid  water  is  transported  out 
from  the  CCL  to  the  cathode  channel.  Similar  to  the  anode  gas 
pressure,  the  cathode  liquid  pressure  is  continuous  across  the 
interfaces  VI  and  VII  while  the  gradients  are  different  due  to 
the  different  diffusivities  (K  and  kr\)  at  different  layers,  which 


Fig.  5.  Distributions  of  (a)  liquid  pressure  in  the  cathode,  and  (b)  liquid  saturation  in  the  cathode  at  different  current  densities. 
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can  be  understood  by  referring  to  Eq.  (17).  With  the  increase 
in  current  density,  the  liquid  pressure  also  increases  due  to  the 
increased  water  generation  from  the  ORR  and  water  crossover 
in  the  CCL.  Since  the  liquid  pressure  in  the  anode  is  nearly  a 
constant  (101.3  kPa),  the  liquid  pressure  in  the  cathode  is  always 
larger,  and  the  induced  back  convection  of  water  from  the  cath¬ 
ode  to  the  anode  is  utilized  to  reduce  the  water  crossover  in  the 
DMFC. 

The  corresponding  distribution  of  cathode  liquid  saturation 
at  different  current  densities  is  shown  in  Fig.  5b.  It  is  seen  that 
cathode  liquid  saturation  is  rather  small  (<0.15)  in  the  cathode 
porous  region,  accordingly  the  mass  transport  resistance  of  oxy¬ 
gen  is  rather  small  (see  Fig.  3b).  The  liquid  saturation  increases 
from  the  cathode  channel  across  the  CDF,  and  similar  to  the 
anode  liquid  saturation,  it  shows  a  discontinuity  at  the  interface 
VII,  which  is  caused  by  the  large  difference  in  permeability  and 
porosity  between  the  CMPF  and  the  CDF.  The  much  smaller 
permeability  of  the  CMPF  mainly  results  in  the  lower  liquid 
saturation  in  this  finer  porous  region.  Similarly,  a  discontinuity 
in  the  liquid  saturation  occurs  at  the  interface  VI,  and  the  less 
hydrophobic  property  ( 0 )  of  the  CCF  leads  to  the  relatively  high 
liquid  saturation.  It  is  also  found  that  the  cathode  liquid  satu¬ 
ration  increases  with  the  increase  in  current  density  due  to  the 
increased  water  flux  across  the  porous  region. 

Comparing  Figs.  4b  and  5b,  it  can  be  seen  that  liquid-gas 
two-phase  flow  takes  place  in  both  the  ACF  and  the  CCF,  and 
thus  the  electrolyte  polymer  is  not  exactly  in  equilibrium  with 
liquid  water  in  the  anode  or  saturated  water  vapor  in  the  cathode. 
Accordingly,  it  is  easy  to  conceive  that  water  content  across 
the  electrolyte  region  is  not  uniform,  and  it  depends  highly  on 
the  liquid  saturation  in  the  ACF  and  the  CCF,  which  in  turn 
influences  the  water  crossover  through  the  membrane. 

The  water  content,  representing  the  dissolved  water  concen¬ 
tration  from  Eq.  (4),  across  the  electrolyte  region  (ACF,  PEM 
and  CCF)  at  different  current  densities  is  shown  in  Fig.  6a. 
Apparently,  the  distribution  of  water  content  across  the  elec¬ 
trolyte  region  is  highly  not  uniform:  the  water  content  near  the 
anode  is  higher  than  that  near  the  cathode.  For  instance,  at  the 
current  density  of  16.0  mA  cm-2,  water  content  of  electrolyte  in 


the  ACF  is  near  to  21,  whereas  in  the  CCF  it  is  only  between 
13  and  15.  This  significant  difference  in  water  content  between 
the  ACF  and  the  CCF  is  mainly  due  to  the  fact  that  liquid  sat¬ 
uration  in  the  ACF  is  relatively  high  (about  0.9  from  Fig.  4b), 
whereas  in  the  CCF  it  is  rather  low  (about  0.12  from  Fig.  5b). 
With  the  increase  in  current  density,  it  is  seen  that  water  con¬ 
tent  decreases  somewhat  in  the  ACF  while  increases  slightly 
in  the  CCF,  leading  to  reduced  gradients  across  the  membrane 
as  well.  This  change  in  water  content  is  also  mainly  caused  by 
the  slight  change  in  liquid  saturation  with  current  density:  liq¬ 
uid  saturation  decreases  in  the  ACF  (Fig.  4b)  while  increases 
in  the  CCF  (Fig.  5b)  with  the  increase  in  current  density.  The 
non-uniform  distribution  of  water  content  and  its  variation  with 
current  density  shown  in  Fig.  6a  clearly  demonstrate  that  water 
content  across  the  electrolyte  region  highly  depends  on  the  liq¬ 
uid  saturation  in  both  the  ACF  and  the  CCF.  It  is  also  noted  that 
near  the  limiting  current  density  (484  mA  cm-2)  water  content 
across  the  ACF  increases  first  from  the  AMPF/ACF  interface 
(III)  and  then  decreases,  which  is  different  from  the  continu¬ 
ous  decrease  at  relatively  low  current  densities.  This  anomalous 
behavior  is  due  to  that  near  the  limiting  current  density,  the  gen¬ 
erated  current  across  the  ACF  is  highly  not  uniform:  current  is 
mainly  generated  near  the  interface  III.  The  generated  current  in 
the  region  near  the  interface  III  can  lead  to  a  considerable  anal¬ 
ogous  dissolved  water  “loss”  by  the  effect  of  electro-osmotic 
drag,  which  can  be  understood  by  referring  to  Eq.  (23). 

Based  on  the  water  content  distribution  shown  in  Fig.  6a, 
we  can  now  examine  how  the  flux  of  water  crossover  varies 
with  current  density  by  studying  Fig.  6b.  It  is  seen  that  the  total 
flux  of  water  crossover  increases  first  slowly  and  then  faster 
with  the  increase  in  current  density,  which  is  consistent  with  the 
experimental  data  reported  elsewhere  [1,2].  As  discussed  in  the 
preceding  section,  the  total  flux  is  due  to  three  water-crossover 
mechanisms:  diffusion  flux,  electro-osmotic  drag  flux  and  back- 
flow  flux.  First,  it  is  seen  that  the  water  flux  by  electro-osmotic 
drag  increases  nearly  linearly  with  current  density.  Secondly, 
the  diffusion  flux  from  the  anode  to  the  cathode  is  quite  sub¬ 
stantial,  which  is  due  to  that  the  water  content  at  the  anode  is 
much  higher  than  that  at  the  cathode  (Fig.  6a).  With  the  increase 
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Distance  across  the  ACL/PEM/CCL  (gm)  Current  density  (mA  cm*2) 


Fig.  6.  Distribution  of  water  content  in  the  electrolyte  region  at  different  current  densities  (a),  and  variation  in  the  water-crossover  flux  through  the  membrane  and 
the  corresponding  net  water- transport  coefficient  with  current  density  (b). 
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in  current  density,  the  diffusion  flux  decreases  as  a  result  of  the 
decreased  water  content  in  the  anode  and  the  increased  water 
content  in  the  cathode,  as  shown  in  Fig.  6a.  Hence,  the  assump¬ 
tion  of  uniform  water  content  across  the  membrane  made  in 
many  previous  papers  on  modeling  DMFCs  will  cause  signifi¬ 
cant  errors  in  predicting  water  crossover  and  cell  performance. 
Third,  the  back-flow  flux  from  the  cathode  to  the  anode,  indi¬ 
cating  by  the  negative  value,  is  slightly  enhanced  with  current 
density,  which  is  due  to  the  increased  liquid  pressure  difference 
across  the  membrane  as  shown  in  Fig.  5a.  Fig.  6b  also  shows 
the  corresponding  net  water-transport  coefficient,  a ,  which  is 
determined  from  Eq.  (39).  It  is  seen  that  a  decreases  rapidly  at 
low  current  densities,  but  the  decrease  becomes  slower  at  high 
current  densities,  from  about  15  to  2.6  when  the  current  density 
is  increased  from  35.0  to  400.0  mA  cm-2.  It  is  worth  mention¬ 
ing  here  that  the  variation  trends  of  the  water-crossover  flux  by 
different  transport  mechanisms,  as  well  as  the  corresponding  a , 
are  all  consistent  with  the  experimental  results  in  our  previous 
work  [1,2]. 


3.2.  Effect  of  contact  angle  of  the  cathode  MPL 


From  the  results  shown  in  Figs.  4-6,  it  is  clear  that  the  prop¬ 
erties  of  cathode  porous  region  have  great  effects  on  water 
crossover  through  the  membrane,  due  to  the  fact  that  the  flux 
of  water  crossover  depends  both  on  the  liquid  pressure  built  up 
in  the  cathode  and  on  the  liquid  saturation  in  the  CCL.  From 
Eqs.  (17)  and  (20),  we  can  obtain: 


Kkr i  a  cos  9(8 K )°'5/trl  d J(s) 

u  = - V  p\  c  = - Vs 

H\/P\  ’  Ml/Pl  ds 


(41) 


which  shows  that  the  liquid  pressure  and  liquid  saturation  in 
the  cathode  can  be  significantly  affected  by  the  permeability, 
contact  angle,  porosity  and  thickness  of  the  cathode  porous  layer, 
especially  the  CMPL.  In  this  section,  the  effect  of  contact  angle 
of  CMPL  is  discussed. 

The  contact  angle  of  the  CMPL  varies  from  110°  to  150°, 
indicating  relatively  low  to  high  hydrophobic  level,  while  all  the 


Table  4 

Limiting  current  densities  under  different  parameters  (rja  =  0.6  V) 


Limiting  current  density  (mAcm  2) 
Effect  of  contact  angle  of  CMPL 


Effect  of  thickness  of  CMPL 


Effect  of  permeability  of  CMPL 


Effect  of  contact  angle  of  CCL 


110° 

484.6 

o 

o 

484.2 

130° 

483.5 

o 

O 

'xl' 

482.6 

150° 

481.9 

15  fxm 

484.5 

30  jxm 

484.2 

60  fxm 

483.6 

90  jxm 

483.1 

5.0e-13m2 

484.7 

2.5-13  m2 

484.2 

1.0e-13  m2 

480.3 

0.75e-13m2 

477.4 

0.5e-13  m2 

470.1 

0.4e-13  m2 

462.8 

0.3e-13  m2 

449.7 

0.25e-13m2 

438.9 

93° 

454.2 

94° 

471.7 

95° 

477.4 

97° 

481.4 

O 

o 

o 

483.2 

o 

o 

483.6 

110° 

483.1 

120° 

481.8 

other  parameters  are  kept  the  same  as  listed  in  Tables  1  and  3. 
Fig.  7a  shows  the  distribution  of  cathode  liquid  pressure  with  dif¬ 
ferent  CMPL  contact  angles  at  nearly  the  same  limiting  current 
densities  (~484.0mAcm-2  as  shown  in  Table  4).  Apparently, 
the  increase  in  the  liquid  pressure  across  the  CMPL  becomes 
more  rapidly  with  the  increase  in  the  contact  angle,  which  also 
leads  to  increased  liquid  pressure  built  up  in  the  CCL.  This  indi¬ 
cates  that  using  highly  hydrophobic  CMPL  helps  build  up  a 
higher  liquid  pressure  in  the  CCL  and  thus  enhance  the  water 
back  flow  from  the  cathode  to  the  anode.  From  Eq.  (41),  the 
higher  liquid  pressure  gradient  in  the  CMPL  with  a  larger  con¬ 
tact  angle  comes  from  the  reduced  liquid  saturation,  as  shown 


Fig.  7.  Distributions  of  (a)  liquid  pressure  in  the  cathode,  and  (b)  liquid  saturation  in  the  cathode  with  different  contact  angles  of  the  CMPL  and  at  the  limiting  current 
densities  (ija  =  0.6  V). 
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Distance  across  the  ACL/PEM/CCL  (jam) 


(b) 


Fig.  8.  Distribution  of  water  content  in  the  electrolyte  region  with  different  contact  angles  of  the  CMPL  and  at  the  limiting  current  densities  (r]a  =  0.6  V)  (a),  and  the 
corresponding  variation  in  the  water-crossover  flux  through  the  membrane  and  the  net  water-transport  coefficient  (b). 


in  Fig.  7b.  It  is  seen  that  the  liquid  saturation  in  the  CMPL 
is  reduced  from  about  0.10  to  0.05  when  the  contact  angle  is 
increased  from  110°  to  150°,  due  to  the  enhanced  hydrophobic 
level.  However,  the  liquid  saturation  in  the  CCL  is  increased 
from  about  0.12  to  0.17,  which  is  caused  by  the  increased  liq¬ 
uid  pressure  in  the  CCL.  This  result  indicates  that  increasing 
the  hydrophobic  level  of  the  CMPL  helps  to  reduce  the  water 
saturation  across  it,  which  is  beneficial  for  the  oxygen  transport 
through  it,  while  the  accordingly  increased  liquid  pressure  built 
up  in  the  CCL  enhances  the  water  flooding  in  the  CL. 

The  water  content  across  the  electrolyte  region  with  differ¬ 
ent  contact  angles  of  the  CMPL  is  shown  in  Fig.  8a.  With  the 
contact  angle  increasing  from  1 10°  to  150°,  water  content  in  the 
CCL  increases  somewhat,  which  is  due  to  the  increased  liquid 
saturation  in  the  CCL,  as  shown  in  Fig.  7b.  As  a  result,  the  gra¬ 
dient  of  water  content  across  the  membrane  is  slightly  reduced. 
Accordingly,  the  change  in  the  cathode  liquid  pressure  and  the 
water  content  in  the  polymer  caused  by  the  change  in  the  CMPL 
contact  angle  will  affect  the  flux  of  water  crossover  from  the 
anode  to  the  cathode,  which  is  shown  in  Fig.  8b.  It  can  be  seen 


that  with  the  increase  in  the  CMPL  contact  angle,  the  water  flux 
by  diffusion  from  the  anode  to  the  cathode  is  slightly  reduced 
due  to  the  lowered  gradient  of  water  content  across  the  mem¬ 
brane,  while  the  water  flux  by  back  flow  from  the  cathode  to 
the  anode  is  slightly  increased  due  to  the  enhanced  liquid  water 
pressure  difference  between  the  cathode  and  the  anode.  As  a 
result,  the  total  flux  of  water  crossover  from  the  anode  to  the 
cathode  as  well  as  the  corresponding  net  water-transport  coeffi¬ 
cient  is  slightly  reduced.  However,  it  is  clear  that  the  reduction 
of  water  crossover  caused  by  the  increased  CMPL  contact  angle 
is  rather  small,  possibly  indicating  that  the  water  crossover  from 
the  anode  to  the  cathode  is  not  very  sensitive  to  the  contact  angle 
of  the  hydrophobic  CMPL. 

3.3.  Effect  of  thickness  of  the  cathode  MPL 

The  effect  of  thickness  of  the  CMPL  is  studied  by  chang¬ 
ing  the  thickness  from  15  to  90  [xm,  while  keeping  all  the 
other  parameters  the  same,  as  listed  in  Tables  1  and  3.  Fig.  9a 
shows  the  distribution  of  cathode  liquid  pressure  with  different 
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Fig.  9.  Distributions  of  (a)  liquid  pressure  in  the  cathode,  and  (b)  liquid  saturation  in  the  cathode  with  different  thicknesses  of  the  CMPL  and  at  the  limiting  current 
densities  (^a  =  0.6  V). 
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Fig.  10.  Distribution  of  water  content  in  the  electrolyte  region  with  different  thicknesses  of  the  CMPL  and  at  the  limiting  current  densities  (^a  =  0.6  V)  (a),  and  the 
corresponding  variation  in  the  water-crossover  flux  through  the  membrane  and  the  net  water- transport  coefficient  (b). 


CMPL  thicknesses  at  nearly  the  same  limiting  current  densities 
(~484.0  mA  cm-2  as  shown  in  Table  4).  It  is  seen  that  the  gra¬ 
dients  of  the  increase  in  the  liquid  pressure  across  the  CMPL 
for  different  CMPL  thicknesses  are  nearly  the  same  due  to  the 
same  properties,  but  a  thicker  CMPL  can  build  up  a  higher  liquid 
pressure  due  to  the  longer  transport  route,  which  also  leads  to  an 
increased  liquid  pressure  built  up  in  the  CCL.  This  indicates  that 
increasing  thickness  of  the  CMPL  is  beneficial  for  building  up  a 
higher  liquid  pressure  in  the  CCL  and  thus  enhancing  the  water 
back  flow  from  the  cathode  to  the  anode.  The  corresponding  liq¬ 
uid  saturation  in  the  cathode  is  shown  in  Fig.  9b.  Similarly,  the 
gradients  of  the  increase  in  the  liquid  saturation  across  the  CMPL 
for  different  CMPL  thicknesses  are  also  the  same,  but  thicker 
CMPLs  can  lead  to  larger  liquid  saturation  near  the  CCL/CMPL 
interface  (VI).  For  instance,  at  the  interface  VI,  liquid  saturation 
in  the  CMPL  for  the  thickness  of  15  p>m  is  about  0.07,  which  is 
smaller  than  that  in  the  CDL  near  the  interface  VII,  while  for 
the  thickness  of  90  jim  it  is  about  0.09,  which  is  larger  than  that 
in  the  CDL  near  the  interface  VII.  It  is  also  found  that  the  liq¬ 
uid  saturation  in  the  CCL  increases  with  the  CMPL  thickness, 
due  to  the  increased  liquid  pressure  in  the  CCL.  The  increased 
liquid  pressure  shown  in  Fig.  9a  and  the  increased  liquid  satura¬ 
tion  in  the  CCL  shown  in  Fig.  9b  will  obviously  affect  the  water 
crossover  through  the  membrane. 

The  water  content  across  the  electrolyte  region  with  differ¬ 
ent  CMPL  thicknesses  is  shown  in  Fig.  10a.  Water  content  in 
the  CCL  is  increased  slightly  due  to  the  increased  liquid  satu¬ 
ration  in  the  CCL,  as  shown  in  Fig.  9b,  which  leads  to  a  slight 
decrease  in  the  gradient  of  water  content  across  the  membrane. 
The  change  in  the  water-crossover  flux  with  the  increase  in  the 
CMPL  thickness  is  shown  in  Fig.  10b.  It  is  also  found  that  the 
water-crossover  flux  by  diffusion  is  slightly  reduced  with  the 
increase  in  the  thickness  due  to  the  lowered  water  content  gradi¬ 
ent  across  the  membrane,  but  the  water-crossover  flux  by  back 
flow  is  slightly  increased  due  to  the  enhanced  liquid  water  pres¬ 
sure  difference  between  the  anode  and  the  cathode.  Accordingly, 
the  total  flux  of  water  crossover  as  well  as  the  corresponding 
net  water-transport  coefficient  is  slightly  reduced.  The  rather 


small  change  in  the  water-crossover  flux  with  the  increase  in 
the  CMPL  thickness  possibly  indicates  that  the  water  crossover 
through  the  membrane  is  also  not  very  sensitive  to  the  thickness 
of  the  hydrophobic  CMPL.  At  the  same  time,  it  should  be  noted 
that  increasing  the  CMPL  thickness  will  lead  to  a  larger  oxygen 
transport  resistance  in  both  the  CMPL  and  the  CCL,  due  to  the 
increased  liquid  saturation,  as  shown  in  Fig.  9b. 

3.4.  Effect  of  permeability  of  the  cathode  MPL 

The  effect  of  permeability  of  the  CMPL  is  examined 
by  changing  its  value  from  5.0  x  10"13  to  0.25  x  1(T13  m2, 
while  keeping  all  the  other  parameters  the  same,  as  listed  in 
Tables  1  and  3.  Fig.  11a  shows  the  distribution  of  cathode 
liquid  pressure  with  different  permeabilities  of  the  CMPL  at 
the  limiting  current  densities  (shown  in  Table  4).  It  is  seen 
that  the  increase  in  the  liquid  pressure  across  the  CMPL  can 
become  much  more  rapidly  with  the  decrease  in  the  CMPL  per¬ 
meability.  For  instance,  the  liquid  pressure  increases  slightly 
from  2.9  to  3.0 kPa  across  the  CMPL  for  the  permeability  of 
5.0  x  10“ 13  m2,  whereas  it  can  increase  about  four  times  from 
about  2.6  to  10.8  kPa  across  the  CMPL  for  the  permeability  of 
0.25  x  10-13  m2.  As  a  result,  the  liquid  pressure  in  the  CCL  is 
significantly  enhanced  with  the  reduction  of  the  CMPL  perme¬ 
ability. 

The  corresponding  liquid  saturation  in  the  cathode  is  shown 
in  Fig.  lib.  It  is  seen  that  with  the  permeability  decreasing 
from  5.0  x  10- 13  to  1.0  x  10  13  m2,  the  liquid  saturation  in  the 
CMPL  is  evidently  decreased.  However,  when  the  permeability 
is  reduced  below  1.0  x  10“ 13  m2,  the  change  in  the  saturation 
in  the  CMPL  is  rather  small.  In  the  CCL,  the  liquid  saturation  is 
increased  significantly  from  about  0. 12  to  0.84  with  the  decrease 
in  the  CMPL  permeability.  This  sharp  increase  in  the  liquid  sat¬ 
uration  is  due  to  the  significantly  enhanced  liquid  pressure  in 
the  CCL,  and  it  obviously  increases  the  risk  of  water  flood¬ 
ing.  The  sharply  increased  liquid  saturation  in  the  CCL  can 
lead  to  two  effects.  On  one  hand,  it  increases  the  oxygen  trans¬ 
port  resistance  through  the  CCL  due  to  the  reduced  gas  void 
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Fig.  11.  Distributions  of  (a)  liquid  pressure  in  the  cathode,  and  (b)  liquid  saturation  in  the  cathode  with  different  permeabilities  of  the  CMPL  and  at  the  limiting 
current  densities  (rja  =  0.6  V). 


fraction,  which  can  be  proved  by  the  reduced  oxygen  concen¬ 
tration  in  the  CCL,  as  shown  in  Fig.  12.  On  the  other  hand,  the 
effective  catalyst  sites  available  for  ORR  reaction  are  seriously 
reduced  due  to  the  increased  coverage  by  liquid  water.  These 
two  effects  result  in  increased  mass  transport  polarization  and 
thus  reduced  cell  performance.  From  Table  4,  it  can  be  seen 
that  the  limiting  current  density  can  be  decreased  from  about 
484.0  to  439.0  mA  cm-2  due  to  the  increased  liquid  saturation  in 
the  CCL. 

The  water  content  across  the  electrolyte  region  with  different 
CMPL  permeabilities  is  shown  in  Fig.  13a.  Due  to  the  highly 
increased  liquid  saturation  in  the  CCL,  it  is  seen  that  water 
content  in  the  CCL  is  also  highly  increased  from  about  13.5 
for  the  permeability  of  5.0  x  10“ 13  m2  to  about  21  for  the  per¬ 
meability  of  0.25  x  10-13m2,  which  is  even  higher  than  that 
in  the  anode  side.  Accordingly,  the  gradient  of  water  content 
across  the  membrane  is  also  reduced  greatly.  The  change  in  the 
water-crossover  flux  with  the  effect  of  the  CMPL  permeability 
is  shown  in  Fig.  13b.  It  can  be  seen  that  the  total  water-crossover 


Distance  across  the  CCL/CMPL/CDL  (pm) 

Fig.  12.  Distributions  of  gas  oxygen  concentration  through  the  cathode  with 
different  permeabilities  of  the  CMPL  and  at  the  limiting  current  densities 
foa  =  0.6  V). 


flux  as  well  as  the  flux  by  diffusion  only  decreases  slightly 
when  the  CMPL  permeability  is  reduced  from  5.0  x  10  13  to 
1.0  x  10-13m2.  However,  when  the  permeability  is  reduced 
below  1.0  x  10-13  m2,  the  water  flux  by  diffusion  is  decreased 
sharply,  and  it  can  even  become  negative,  indicating  diffusing 
from  the  cathode  to  the  anode.  This  rapidly  change  in  the  diffu¬ 
sion  flux  is  caused  by  the  significantly  increased  water  content 
in  the  CCL.  On  the  other  hand,  the  highly  enhanced  liquid 
pressure  in  the  CCL  leads  to  a  rapid  increase  in  the  flux  of 
back  flow  from  the  cathode  to  the  anode.  As  a  result,  when  the 
permeability  is  reduced  from  1.0  x  10“13  to  0.25  x  10“13m2, 
the  total  water-crossover  flux  can  decrease  greatly  from 
about  11.5  to  3.3  juimolcm-2  s-1,  and  the  corresponding  net 
water-transport  coefficient  is  reduced  significantly  from  about 
2.4  to  0.8. 

Apparently,  decreasing  the  permeability  of  the  CMPL  is  a 
very  effective  way  to  build  up  high  liquid  pressure  in  the  CCL 
and  thus  reduce  the  water  crossover  from  the  anode  to  the  cath¬ 
ode.  However,  the  induced  relatively  high  liquid  saturation  in  the 
CCL  is  not  beneficial  for  the  oxygen  transport  and  cell  perfor¬ 
mance.  In  fact,  the  liquid  saturation  in  the  CCL  can  be  adjusted 
by  changing  the  hydrophobic  level  or  the  contact  angle  of  the 
CCL,  as  discussed  in  the  next  section. 

3.5.  Effect  of  contact  angle  of  the  cathode  CL 

This  section  presents  the  effect  of  contact  angle  of  the  CCL  on 
water  crossover.  The  contact  angle  is  changed  from  93°  to  120°, 
while  the  permeability  of  the  CMPL  is  fixed  at  0.75  x  l(r13m2 
and  all  the  other  parameters  are  kept  the  same  as  listed  in 
Tables  1  and  3.  Fig.  14a  shows  the  distribution  of  cathode  liquid 
pressure  with  different  contact  angles  of  the  CCL  at  the  limiting 
current  densities  (shown  in  Table  4).  It  is  seen  that  the  increase 
in  the  liquid  pressure  across  the  CCL  becomes  more  rapidly 
with  the  increase  in  the  contact  angle  of  the  CCL.  For  instance, 
the  liquid  pressure  increases  negligibly  from  5.70  to  5.71  kPa 
across  the  CCL  for  the  contact  angle  of  93°,  but  it  can  increase 
over  2  times  from  about  6.0  to  13.5  kPa  across  the  CCL  for  the 
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Fig.  13.  Distribution  of  water  content  in  the  electrolyte  region  with  different  permeabilities  of  the  CMPL  and  at  the  limiting  current  densities  (^a  =  0.6  V)  (a),  and 
the  corresponding  variation  in  the  water-crossover  flux  through  the  membrane  and  the  net  water-transport  coefficient  (b). 


contact  angle  of  120°.  This  indicates  that  the  liquid  pressure  in 
the  CCL  can  also  be  significantly  enhanced  by  increasing  the 
hydrophobic  level  of  the  CCL. 

The  corresponding  liquid  saturation  in  the  cathode  is  shown 
in  Fig.  14b.  With  the  increase  in  the  contact  angle  from  93°  to 
110°,  the  liquid  saturation  in  the  CCL  is  significantly  decreased, 
but  the  decrease  in  the  liquid  saturation  becomes  rather  slow 
when  the  contact  angle  is  above  110°.  As  such,  increasing  the 
hydrophobic  level  of  the  CCL  can  not  only  effectively  build 
up  high  liquid  pressure,  which  is  beneficial  for  reducing  water 
crossover,  but  also  helps  reduce  the  water  flooding  in  the  CCL. 
It  can  be  seen  from  Table  4  that  the  limiting  current  density  can 
be  increased  from  about  454.0  to  482.0  mA  cm-2  resulting  from 
the  reduced  liquid  saturation. 

The  water  content  across  the  electrolyte  region  with  different 
contact  angles  of  the  CCL  is  shown  in  Fig.  15a.  Due  to  the 
highly  decreased  liquid  saturation  in  the  CCL,  water  content  in 
the  CCL  is  also  significantly  decreased  from  about  19.3  for  the 


contact  angle  of  93°  to  about  13.5  for  the  contact  angle  of  120°. 
Accordingly,  the  gradient  of  water  content  across  the  membrane 
is  also  increased  greatly.  The  increased  liquid  pressure  and  the 
decreased  liquid  saturation  with  the  contact  angle  in  the  CCL 
result  in  an  inconsistent  effect  on  the  flux  of  water  crossover, 
which  is  shown  in  Fig.  15b.  It  is  interesting  to  see  that  increasing 
the  contact  angle  of  the  CCL  can  simultaneously  enhance  the 
back- flow  flux  from  the  cathode  to  the  anode,  which  is  beneficial 
for  reducing  water  crossover,  but  increase  the  diffusion  flux  from 
the  anode  to  the  cathode,  which  tends  to  increase  water  crossover. 
This  inconsistent  effect  leads  to  the  result  that  an  intermediate 
contact  angle  of  about  97°  yields  the  largest  total  flux  of  water 
crossover,  and  decreasing  or  increasing  the  contact  angle  can  all 
effectively  reducing  the  total  flux  of  water  crossover.  However, 
from  the  point  view  of  oxygen  transport  and  cell  performance, 
increasing  the  contact  angle  of  the  CCL  is  preferred  for  both 
lowered  water  crossover  through  the  membrane  and  improved 
cell  performance. 


Distance  across  the  CCL/CMPL/CDL  (pm) 


Fig.  14.  Distributions  of  (a)  liquid  pressure  in  the  cathode,  and  (b)  liquid  saturation  in  the  cathode  with  different  contact  angles  of  the  CCL  and  at  the  limiting  current 
densities  (ija  =  0.6  V). 
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Fig.  15.  Distribution  of  water  content  in  the  electrolyte  region  with  different  contact  angles  of  the  CCL  and  at  the  limiting  current  densities  (^a  =  0.6  V)  (a),  and  the 
corresponding  variation  in  the  water-crossover  flux  through  the  membrane  and  the  net  water- transport  coefficient  (b). 


4.  Conclusions 

In  this  work,  we  present  a  one-dimensional,  isothermal 
two-phase  mass  transport  model  to  investigate  the  water  trans¬ 
port  through  the  membrane  electrode  assembly  for  liquid-feed 
DMFCs.  Liquid-gas  two-phase  mass  transport  in  the  porous 
anode  and  cathode  is  modeled  with  classical  multiphase  flow 
theory  in  porous  media,  in  which  the  liquid  phase  consists  of 
methanol-water  solution,  while  the  gas  phase  includes  carbon 
dioxide  gas,  methanol  vapor,  and  water  vapor.  In  the  electrolyte 
region  that  includes  the  membrane  and  the  catalyst  layers  (CLs), 
transport  of  dissolved  water  is  driven  by  diffusion,  electro- 
osmotic  drag  and  convection.  In  particular,  water  transport  in 
three  phases  (i.e.,  liquid,  vapor  and  dissolved  phase)  in  the 
CLs  is  coupled  intrinsically,  and  finite-rate  interfacial  exchanges 
between  different  phases  are  introduced  without  any  phase  equi¬ 
librium  assumptions.  With  this  new  model,  we  investigated 
each  water-crossover  mechanism  through  the  membrane  for  the 
DMFC,  and  studied  the  effect  of  cathode  structure,  such  as  the 
design  of  the  cathode  MPL,  on  the  water  crossover  through 
the  membrane  and  the  water  ejection  through  the  cathode.  The 
salient  findings  in  this  work  are  summarized  as  follows: 

( 1 )  At  low  current  densities,  liquid  saturation  in  the  anode  CL  is 
in  general  larger  than  that  in  the  cathode,  resulting  in  a  large 
difference  in  water  content  across  the  membrane  and  thus 
water  diffusion  from  the  anode  to  the  cathode.  However,  on 
the  other  hand,  liquid  pressure  in  the  cathode  is  larger  than 
that  in  the  anode,  leading  to  back  convection  of  liquid  water 
directly  from  the  cathode  to  the  anode.  Lowering  water  dif¬ 
fusion  and  enhancing  water  back  convection  through  the 
membrane  should  be  utilized  to  reduce  the  water  crossover 
in  the  DMFC. 

(2)  With  the  increase  in  current  density,  liquid  saturation 
decreases  in  the  anode  CL,  but  increases  in  the  cath¬ 
ode  CL.  At  the  same  time,  an  increase  in  current  density 
also  leads  to  an  increase  in  liquid  pressure  in  the  cath¬ 
ode  CL.  Accordingly,  the  water-crossover  flux  by  diffusion 


decreases  significantly,  although  the  flux  by  back  convection 
increases  to  some  extent.  As  a  result,  the  net  water-transport 
coefficient  decreases  rapidly  at  low  current  densities,  but  the 
decrease  becomes  slower  at  high  current  densities. 

(3)  It  is  found  that  increasing  the  liquid  pressure  in  the  cathode 
CL  by  lowering  the  permeability  or  increasing  the  contact 
angle  or  the  thickness  of  the  cathode  MPL  is  accompanied 
with  an  increase  in  the  liquid  saturation  in  it,  which  helps 
reduce  the  water  diffusion  from  the  anode  to  the  cathode. 
As  a  result,  water  crossover  can  be  reduced  significantly 
by  lowering  the  permeability  of  the  cathode  MPL  due  to 
the  enhanced  back  convection  and  decreased  water  diffu¬ 
sion,  while  the  reduction  of  water  crossover  by  increasing 
the  hydrophobic  level  or  the  thickness  of  the  cathode  MPL 
is  relatively  small.  Morever,  lowering  the  permeability  or 
increasing  the  hydrophobic  level  of  the  cathode  MPL  can 
reduce  the  water  saturation  in  it,  which  is  beneficial  for 
the  oxygen  transport.  However,  it  should  be  noted  that  the 
induced  relatively  high  liquid  saturation  in  the  cathode  CL 
increases  the  risk  of  electrode  flooding.  The  liquid  satura¬ 
tion  in  the  cathode  CL  can  be  lowered  by  increasing  its 
hydrophobic  level,  which  is  also  very  effective  to  build 
up  high  liquid  pressure  and  thus  reduce  the  flux  of  water 
crossover. 
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